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Abstract 

An exact analytical diagonalization is used to solve the two dimensional Ex- 
tended Hubbard Model for system with finite size. We have considered an Extended 
Hubbard Model (EHM) including on-site and off-site interactions with interaction 
energy U and V respectively, for square lattice containing 4x4 sites at one-eighth 
filling with periodic boundary conditions, recently treated by Kovacs et al pp. 
Taking into account the symmetry properties of this square lattice and using a 
translation linear operator, we have constructed a r — space basis, only, with 85 
state- vectors which describe all possible distributions for four electrons in the 4x4 
square lattice. The diagonalization of the 85 x 85 matrix energy allows us to study 
the local properties of the above system as function of the on-site and off-site 
interactions energies, where, we have shown that the off-site interaction encour- 
ages the existence of the double occupancies at the first exited state and induces 
supplementary conductivity of the system. 
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1 Introduction 



In the last years, the Extended Hubbard Model (EHM) was introduced to explain some 
interesting phenomena including metal-insulator transition [2J [31 HI EJ [6] , antiferromag- 
netism [3 El El ED] , and high-T c superconductivity pH H21 [Lj . This EHM is a standard 
simple model of interacting itinerant electrons in a solid [HI [151 EE El HH1 ?]■ Although 
this model is too idealized to be regarded as a quantitatively reliable model of real solids, 
it contains physically essential features of interacting itinerant electron systems. Despite 
being one of the most studied models to describe strongly correlated electrons system, 
many questions concerning the Hubbard model remain as open problems. In [19J, we 
have applied the Self Consistent Random Phase Approximation (SCRPA) [201 EI] to 
solve the EHM in ID, where, we have shown that this approach treat the correlations 
of closed chains with a rigours manner. The behaviour of our SCRPA ground state 
and gap energies shows that the repulsive off-site interaction between the electrons of 
the neighbouring atoms induces supplementary conductivity, since, the SCRPA energy 
gap vanishes when the closed chains of the EHM are governed by a strong repulsive on- 
site and an intermediate repulsive off-site interactions. But, due to the restricted motion 
along one direction in space, the Hubbard model chain doesn't exhibit any ferromagnetic 
feature. Thus, the ID Hubbard is a nice prototype to describe, only, the ID classical 
conductors. 

Stimulated by the discovery of high-T c superconductivity in the cuprate plans, the 
2D Hubbard model have attracted a great attention in the recent two decades [221 123] . In 
spite of its simple description in square lattice, it is not obvious to solve the 2D Hubbard 
model in the general case. The exact solution of 2D Hubbard model is not still reached, 
but, a great variety of approximate treatments have been proposed [2H 125] . 

Very recently, E. Kovacs et al. [1] proposed an exact solution of the Usual Hubbard 
Model for 4sites x Asites cluster with low concentration, where each cluster presents four 
electrons. The dynamics of these electrons is described by the usual kinetic of electrons 
with hopping energy t and the repulsive on-site interaction between electrons in the same 
site with interaction energy U. 

In the present paper, we purpose to apply this method to an extended Hubbard model 
which takes into account the off-site interaction with an interaction energy V 

This paper is organized as follows. In sec. 2, we present the model and the calculation 
procedure which allows us to construct the Hamiltonian matrix of dimension 85 x 85. In 
sec. 3, we present our results for ground state and gap energies and discuss the effect of the 
off-site interaction on the dynamics of system. Finally, in Sec. 4 we give our conclusions. 



2 



2 Model and Formalism 



The two-dimensional Extended Hubbard Model on a square lattice is given by: 

H = - 1 ,.!<■]. ,,<■;,. o + U ^n i)t ni4 + V ^ n ij0 -n Ji0 -/ (1) 

(ij>,CT i (i,j),a,a' 

where c- CT (cj jCr ) are the creation (annihilation) operators for a fermion of spin a at 
site i ( j ) with periodic boundary conditions. Thus, t is the hopping term from the site 
j to the site i, where sums over nearest neighbour sites. The second term describes 
the local repulsive interaction with parameter U. The last term takes into account the 
nearest-neighbour repulsion between electrons with energy V. 

The resolution of the model (1) in case of finite size system, gives the exact solution of 
some physical quantities as the ground state energy, the energy gap and the occupation 
numbers. We consider, thus, a two dimensional L x L = N = 16 square lattice at 
one-eighth filling (four electrons per cluster), with periodic boundary conditions in both 
directions. For this considered system, three types of particle configurations may occur. 
First, we can have two double occupancies at sites i and j (i ^ j). Second, we may have 
a double occupancy at site i and two electrons with opposite spins at sites j and k (i ^ j, 
i k and j ^ k). Finally, we may have four single occupancies placed on different sites 
of this 4x4 square lattice. These three possible configurations provide, respectively, the 
three states: 

i a > = (%A-°) ( c J> c i,-*) i°) 

l C ) = ( C *> C i,a) { C l,-a C l-a) 1°) 

where |0) represents the state vacuum with no electron present. 

The states \a) , |6) and |c) generate N d = 16xl !^ 5x15 = 28800 states, which describe 
all possible distributions of our four electrons in the 4x4 cluster. In order to construct a 
r — space basis, it is convient to define a linear operator T pQ , which verifies the relation: 

T(A + B) =T(A)+T(B) 

Where A and B represent the particle configurations as mentioned in Fig. 1, T(A) 
is the linear combination of the 16 contributions obtained by the translation of the 
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Figure 1: The three possible types of particle configurations (a: two double occupancies 
- b: double occupancy and two single occupancies - c: four single occupancies) 

configuration A to each site of the 4x4 cluster. Using this linear operator and taking 
into account the symmetry proprieties of the 4x4 cluster in r — space representation, we 

can regroup these states in 85 cluster states denoted by \n), where n = 1,2, ,85, 

and are all orthogonal vectors. For example, |1) is the linear combination of all states 
| a) with \Ri — Rj\ = a, where Ri (Rj) is the lattice position of the site i (j) and a is the 
square lattice parameter. |2) is the linear combination of all states \b) with \R\ — Rj\ = 
\Rj — Rk\ = a and \Ri — Rk\ = 2a. Whereas |3) is ,also, the linear combination of all 
states 1 6) with \Ri — Rj\ = \Rj — R^\ = a but — R^\ = \/2a (For the definition of 
the other 82 vectors, see ref. PQ). 

The application of the Hamiltonian (1) on the basis vector |1) gives: 

H\l) = (2U + 4V)\l)+t(\2)-\3)) 

It is clear that the Hamiltonian would not be diagonal in our r — space basis, since 
the application of the kinetic term H = — Yli^ij c \ a c j,o on a vector |n) gives, always, the 
new states \n') after the creation and the annihilation of electrons in the different lattice 
sites. Thus, it is necessary to define the matrix energy E 85 x 85 as: 

E nm = (n \H\m) 

Where \n) and |m) are two vectors of the r — space basis. In order to study the local 
properties of the 4x4 square lattice, it is appropriate to determine , numerically, the 
eigenvalues and the eigenvectors of the matrix energy E 85x85 . We consider the obtained 
first and second minimums of the eigenvalues as, respectively, the ground state and the 
first excited energies. 
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3 Results and discussion 



First, we disregarded the off-site interaction, and we have plotted in Fig.2 and Fig. 3 
the ground state and the first excited state energies, respectively, as function of the on- 
site interaction energy U. The corresponding curves of this case (V = 0) show that the 
ground state energy has smooth (less than linear) U dependence, whereas the first excited 
state energy is U independent and fixed at — 8t for any value of U. Thus, it is clear that 
the first excited state is an eigenvector of the kinetic energy <7*ij c l,o- c .?>) > si nce 

this excited state avoid totally the double occupancy Whereas, the ground state is not 
eigenstate of the kinetic energy or the on-site interacting part of H only, but is eigenvector 
for the sum of both. Thus, at the ground state, the dynamics of the electrons system is 
governed by a competition between the habitual kinetic and the on-site interaction. 

Then, we have taken into account the off-site interaction (V ^ 0) and we have plotted 
in in the same previous figures (Fig.2 and Fig. 3) the ground state and the first excited 
state energies, respectively, as function of the on-site interaction energy U for two values 
of V/t . 



Figure 2: Ground state energy as function of U/t for different values of V/t 

The Fig. 2 shows that the energy E has a smooth U dependence and a linear V 
dependence. This behaviour is similar to the one obtained for the ground state energy of 
the chains Extended Hubbard Model with the Self Consistent Random Phase Approxi- 
mation (SCRPA) [26], where we have also, defined the matrix energy but in the basis of 
vectors impulsion-space (k — space) and not in r — space representation. The curves of 
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Figure 3: First excited state energy as function of U/t for different values of V/t 



Fig. 2 show, also, that E decreases with V for a fixed value of U. Thus, we can conclude 
that this off-site interaction imposes the electron system to avoid partially the double 
occupancy in this ground state. 

The Fig. 3 shows that the first exited state energy becomes U dependent for V ^ 0. 
Thus, the off-site interaction encourages the existence of the double occupancies in this 
exited state. For a weak off-site interaction (V/t = 0.5), the corresponding curve shows 
that E* becomes U independent for the high values of U. But for an intermediate off-site 
interaction (V/t = 1), E* still remains dependent U; since, we have the opportunity to 
have the double occupancies even for high values of U. 

In order to analyze the on-site and off-site interactions effects on the repartition of 
our four electrons in the above system, we define the double occupancy coefficient in the 
first exited state D* as the probability to have a couple of electrons (tl) on the same 
site in this exited state, 

D * = ^ ( n *r n a) 

i 

where the sum over all 4 x 4 cluster sites and the mean values (...) are taken in the 
corresponding eigenvector to E*. 

The double occupancy coefficient D* in the first exited state is shown in Fig.4 as 
function of U/t for different values of V/t. For V/t = 0, we have D* = 0. Thus, 
effectively, our system avoids completely the double occupancy at this exited state. But, 
after taking into account the off-site interaction, we have D* ^ 0. For V/t = 1, we have 
D* 7^ even for the strong values of U/t. But, for V/t = 0.5, the coefficient D* vanishes 
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Figure 4: Double occupancy coefficient D* in the first exited state as function of U/t for 
different values of V/t 

for the strong value of U/t, where the off-site interaction becomes very weak before the 
on-site interaction. Thus, the behaviour of our system in this regime (U <C V) is similar 
to the one found in Ref. [1] . Where the authors have shown that 40% of the excited 
states of the 4x4 cluster are U independent. But, it is clear that the number of this U 
independent exited states decreases if we take into account the off-site interaction, since 
this interaction encourage the formation of the double occupancies in the exited states. 

Finally, in order to analyze the effect of the off-site interaction on the dynamics of 
electrons, we define the energy gap Ae as the difference between the first excited state 
energy E* and the ground state energy E 

Ae = E* — E 

In Fig. 5, we plot the variation of this energy gap Ae as function of the repulsive 
on-site interaction energy U for different values of the off-site interaction energy V. 

For a fixed value of V/t, the energy gap decreases with U. We deduce that the 
repulsive on-site interaction (U y 0) increases the conductivity of the system, since the 
repulsion between the two electrons of the same site encourages every electron to jump to 
the neighbouring site. Moreover, this curves show that the off- site interaction increases 
also the conductivity of this 4x4 cluster. For a weak off-site interaction (V/t ^ 0.5), 
the effect of V is remarkable only for the weak on-site interaction whereas it becomes 
practically non-existent for strong on-site interaction. But, for an intermediate off-site 
interaction (V/t = 1), the effect remains remarkable even for the strong values of U/t. 
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Figure 5: Energy gap as function oiU/t for different values of V/t 



Thus, we can conclude that with a strong on-site and an intermediate off-site interactions, 
our 4x4 cluster have an important conductivity. 

4 Conclusion 

In this paper, an exact diagonalization in the r — space was proposed to solve the two 
dimensional Extended Hubbard Model for finite size system. In particular, we have 
considered a two dimensional L x L = N = 16 square lattice at one-eighth filling with 
respecting the periodic boundary conditions in both directions. 

First, the numeric diagonalization of our Hamiltonian matrix allow us to determine 
some interesting local properties of our 4x4 square lattice as: the ground state energy 
Eqs, the first excited state energy E*, the gap energie Ae and the double occupation 
number per site D*. Then, the analysis of the behaviours of these obtained local prop- 
erties as function of U and V allows as to study the distribution and the dynamics of 
the elctrons system in two interesting states: ground and first exited states. Thus, we 
have found that our system has always the double occupancies in the ground state for 
any value of U and V. In the first exited state, we have shown that the off-site site 
interaction encourages the electrons system to form the double occupancies, where the 
coefficient D* vanishes for V = 0. Whereas, for an intermediate off-site interaction, we 
found that we have always the probability to have these double occupancies. But, this 
probability vanishes for strong on-site and weak off-site interactions, where, the electrons 
system avoids completely the double occupancy. Thus, the behaviour of our system in 
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this regime (V <C U) is similar to the one found in Ref. pQ. Finally, the analysis of the 
off-site interaction effect on the energy gap shows that the repulsive off-site interaction 
induces supplementary conductivity of the system, where, this effect of V is more re- 
markable for an intermediate off-site interaction, since we have a reduction of order 20% 
in this regime. 
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